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^ We review recent results obtained in the physics of the thermal Casimir force acting 

between two dielectrics, dielectric and metal, and between metal and semiconductor. 
The detailed derivation for the low-temperature behavior of the Casimir free energy, 
pressure and entropy in the configuration of two real dielectric plates is presented. For 
dielectrics with finite static dielectric permittivity it is shown that the Nernst heat the- 
orem is satisfied. Hence, the Lifshitz theory of the van der Waals and Casimir forces is 
demonstrated to be consistent with thermodynamics. The nonzero dc conductivity of 
dielectric plates is proved to lead to a violation of the Nernst heat theorem and, thus, 
is not related to the phenomenon of dispersion forces. The low-temperature asymptotics 
of the Casimir free energy, pressure and entropy are derived also in the configuration of 
one metal and one dielectric plate. The results are shown to be consistent with thermo- 
C "3 ' dynamics if the dielectric plate possesses a finite static dielectric permittivity. If the dc 

i conductivity of a dielectric plate is taken into account this results in the violation of the 

, Nernst heat theorem. We discuss both the experimental and theoretical results related to 

the Casimir interaction between metal and semiconductor with different charge carrier 
density. Discussions in the literature on the possible influence of spatial dispersion on 
" the thermal Casimir force are analyzed. In conclusion, the conventional Lifshitz theory 

, taking into account only the frequency dispersion remains the reliable foundation for the 

■ interpretation of all present experiments. 
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1. INTRODUCTION 

The Casimir effectP is the force and also the specific polarization of the vacuum 
arising in restricted quantization volumes and originating from the zero-point oscil- 
lations of quantized fields. This force acts between two closely spaced macrobodies, 
between an atom or a molecule and macrobody or between two atoms or molecules. 
During more than fifty years, passed after the discovery of the Casimir effect, it has 
attracted much theoretical attention because of numerous applications in quantum 
field theory, atomic physics, condensed matter physics, gravitation and cosmology, 
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mathematical physics, and in nanotechnology (see monographs [2]-[4] and reviews 
[5HZ1)- In multidimensional Kaluza-Klein supergravity the Casimir effect was usea® 
as a mechanism for spontaneous compactification of extra spatial dimensions a nd 
for constraining the Yukawa-type corrections to Newtonian gravity. ^ 1 1 1 1 1 1 2 J13| j n 
quantum chromodynamics the Casimir energy plays an important role in the bag 
model of hadrons. 1 ^ 1 In cavity quantum electrodynamics the Casimir interaction be- 
tween an isolated atom and a cavity wall leads to the level shifts of atomic electrons 
depending on the position of the atom near the wall!* Both the van der Waals and 
Casimir forces are usecPSGSl f or t ne theoretical interpretation of recent experiments 
on quantum reflection and Bose-Einstein condensation of ultracold atoms on or near 
the cavity wall of different nature. In condensed matter physics the Casimir effect 
turned out to be important for interaction of thin films, in wetting processes, and in 
the theory of colloids and lattice defects!^ The Casimir force was used to actuate 
nanoelectromechanic al d eviceJ^ and to study the absorption of hydrogen atoms 
by carbon nanotubesP^ Theoretical work on the calculation of the Casimir ener- 
gies and forces stimulated important achievements in mathematical physics and in 
the theory of renormalizations conne cted with the method of generalized zeta func- 
tion and heat kernel expansion! ® ! ^ ! All this made the Casimir effect the subject of 
general interdisciplinary interest and attracted permanently much attention in the 
scientific literature. 

The last ten years were marked by the intensive experimental investigation of 
the Casimir force between metallic test bodies (see Refs. |2T)H3"3"|) . During this time 
the agreement between experiment and theory on the level of 1-2% of the measured 
force was achieved. This has become possible due to the use of modern laboratory 
techniques, in particular, of atomic force microscopes and micromechanical torsional 
oscillators. Metallic test bodies provide advantage in comparison with dielectrics be- 
cause their surfaces avoid charging. In Refs. [34j [35l where the importance of the 
Casimir effect for nanotechnology was pioneered, it was demonstrated that at sepa- 
rations below 100 nm the Casimir force becomes larger than the typical electrostatic 
forces acting between the elements of microelectromechanical systems. Bearing in 
mind that the miniaturization is the main tendency in modern technology, it be- 
comes clear that the creation of new generation of nanotechnological devices with 
further decreased elements and separations between them would become impossible 
without careful account and calculation of the Casimir force. 

Successful developments of nanotechnologies based on the Casimir effect calls 
for more sophisticated calculation methods of the Casimir forces. Most of theoreti- 
cal output produced during the first decades after Casimir 's discovery did not take 
into account experimental conditions and real material properties of the boundary 
bodies, such as surface roughness, finite conductivity and nonzero temperature. The 
basic theory giving the unified descriptio n of both the van der Waals and Casimir 
forces was elaborated by 

Lifshit j36l37l3a 

shortly after the publication of Casimir's 
paper. It describes the boundary bodies in terms of the frequency dependent dielec- 
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trie permittivity s(lo) at nonzero temperature T. In the applications of the Lifshitz 
theory to dielectrics it was supposed that the static dielectric permittivity (i.e., the 
dielectric permittivity at zero frequency) is finite. The case of ideal metals was ob- 
tained from the Lifshitz theory by using the so-called Schwinger's prescription)^ 
i.e., that the limit \s(u>)\ — > oo should be taken first and the static limit w — > 
second. For ideal metals the same result, as follows from the Lifshitz theory com- 
bined with the Schwinger's prescription, was obtained independently in the frame- 
work of thermal quantum field theory in Matsubara formulation!^ However, the 
cases of real dielectrics and metals (which possess some nonzero dc conductivity 
at T > and finite dielectric permittivity at nonzero frequencies, respectively) re- 
mained practically unexplored for a long time. The case of semiconductor boundary 
bodies was also unexplored despite of the crucial role of semiconductor materials in 
nanotechnology. 

Starting in 2000, several theoretical groups in different countries attempted to 
describe the Casimir interaction between real metals at nonzero temperature in the 
framework of Lifshitz theory. They have used different models of the metal con- 
ductivity and arrived to controversial conclusions. In Ref. [¥TJ using the dielectric 
function of the Drude model, quite different results than for ideal metals were ob- 
tained. According to Ref. SJ at short separations (low temperatures) the thermal 
correction to the Casimir force acting between real metals is several hundred times 
larger than between ideal ones. In addition, at large separations of a few microm- 
eters (high temperatures) a two times smaller magnitude of the thermal Casimir 
force was fo und than between ideal metals (the latter is known as "the classical 
limit"S^E31) . In Refs. [44] [45] the dielectric permittivity of the plasma model was 
used to describe real metals and quite different results were obtained. At short sep- 
arations the thermal correction appeared to be small in qualitative agreement with 
the case of ideal metals. At large separations for real metals the familiar classical 
limit was reproduced. Later the approach of Ref. [41] was supported in Refs. [46] 
l47l The plasma model approach can be used at such separations a that the char- 
acteristic frequency c/(2a) belongs to the region of infrared optics . Late r a more 
general framework, namely the impedance approach was sugg 

estec™ 

which is 

applicable at any separation larger than the plasma wave length. It was supported 
in Refs. 7, 50, 51] In the region of the infrared optics, the impedance approach leads 
to practically the same results as the plasma model approach. As was shown in 
Refs. [52] 153] the Drude model approach leads to the violation of the Nernst heat 
theorem when applied to perfect metal crystal lattices with no impurities. This ap- 
proach was also excluded by experiment at 99% confidence in the separation region 
from 300 to 500 nm and at 95% confidence in the wider separation region from 170 
to 700 nm P^^ B On the contrary, the plasma model and impedance approaches 
were shown to be in agreement with thermodynamics and consistent with exper- 
iment. The polemic between different theoretical approaches to the description of 
the thermal Casimir force in the case of real metals can be found in Refs. 147 ] I54H561 
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These findings on the application of the Lifshitz theory to real metals have in- 
spired a renewed interest in the Casimir force between dielectrics. As was mentioned 
above, at nonzero temperature dielectrics possess an although small but not equal 
to zero dc conductivity. In Ref. [57] [58] the van der Waals force arising from the dc 
conductivity of a dielectric plate was shown to lead to large effect in noncontact 
atomic frictionJ^Sl a phenomenon having so far no satisfactory theoretical explana- 
tion. This brings up the question: Is it necessary or possible to take into account 
the dc conductivity of dielectrics in the Lifshitz theory? Recall that in the case of 
a positive answer the static dielectric permittivity of a dielectric material would be 
infinitely large. It is amply clear that the resolution of the above issue should be 
in accordance with the fundamentals of thermodynamics. For this reason, it is de- 
sirable to investigate the low-temperature behavior of the Casimir free energy and 
entropy for two dielectric plates both with neglected and included effects of the dc 
conductivity. 

A major breakthrough in the investigation of this problem was achieved in the 
year 2005. In Ref. [SD]a new variant of perturbation theory was developed in a small 
parameter proportional to the product of the separation distance between the plates 
and the temperature. As a result, the behavior of the Casimir free energy, entropy 
and pressure at low temperatures was found analytically. If the static dielectric 
permittivity is finite, the thermal correction was demonstrated to be in accordance 
with thermodynamics. This solves positively the fundamental problem about the 
agreement between the Lifshitz theory and thermodynamics for the case of two 
dielectric plates. In Ref. [60] it was shown that, on the contrary, the formal inclusion 
of a small conductivity of dielectric plates at low frequencies into the model of 
their dielectric response leads to a violation of the Nernst heat theorem. This result 
gives an important guidance on how to extrapolate the tabulated optical data for 
the complex refractive index to low frequencies in numerous applications of the 
van der Waals and Casimir forces. All these problems and related ones arising for 
semiconductor materials are discussed in this review. 

In Sec. 2 we derive the analytical behavior of the Casimir free energy, entropy and 
pressure in the configuration of two parallel dielectric plates at both low and high 
temperature. It is demonstrated that if the static dielectric permittivity is finite the 
Lifshitz theory is in agreement with thermodynamics. Sec. 3 contains the derivation 
of the low-temperature behavior for the Casimir free energy and entropy between 
two dielectric plates with included dc conductivity. In this case the Lifshitz theory 
is found to be in contradiction with the Nernst heat theorem. The conclusion is 
made that the conductivity properties of a dielectric material at a constant current 
are unrelated to the van der Waals and Casimir forces and must not be included 
into the model of dielectric response. In Sec. 4 we consider the thermal Casimir 
force between metal and dielectric. This problem was first investigated in Ref. 1611 
It was found that for dielectric plate with finite static dielectric permittivity the 
Nernst heat theorem is satisfied but the Casimir entropy may take negative values. 
Here we not only reproduce an analytical proof of the Nernst heat theorem but also 
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find the next perturbation orders in the expansion of the Casimir free energy and 
entropy in powers of a small parameter. The results obtained in high-temperature 
limit are also provided. Sec. 5 is devoted to the Casimir interaction between metal 
and dielectric plates with included dc conductivity of the dielectric material. We 
demonstrate that in this case the Nernst heat theorem is violated. Sec. 6 contains 
the discussion of semiconductors which present a wide variety of electric properties 
varying from metallic to dielectric. We consider the Casimir interaction between 
metal and semiconductor test bodies and formulate the criterion when it is appro- 
priate to include the dc conductivity of a semiconductor into t he model of dielectric 
response. In doing so, the results of recent experiment j62 | 63 J on ^ e measurement 
of the Casimir force between metal and semiconductor test bodies are taken into 
account. To this point the assumption has been made that metal, dielectric and 
semiconductor materials of the Casimir plates possess only temporal dispersion, 
i.e., can be described by the dielectric permittivity depending only on frequency. In 
Sec. 7 we discuss recent controversial results by different authors (see, e.g., Refs. [Ml 
- [55]) attempting to take into account also spatial dispersion. As is shown in this 
section, the way of inclusion of spatial dispersion into the Lifshitz theory, used in 
Refs.[Ml-IB51 is unjustified. We argue that the account of spatial dispersion cannot 
influence theoretical results obtained with the help of usual, spatially local, Lifshitz 
theory within presently used ranges of experimental separations. Sec. 8 contains our 
conclusions and discussion. 



2. NEW ANALYTICAL RESULTS FOR THE THERMAL 
CASIMIR FORCE BETWEEN DIELECTRICS 

We consider two thick dielectric plates (semispaces) described by the frequency- 
dependent dielectric permittivity e{u)) and restricted by the parallel planes z — ±a/2 
with a separation a between them, in thermal equilibrium at temperature T. The 
Lifshitz formula for the free energy of the van der Waals and Casimir interaction 
between the plates is given b y ! 2 l 3 l 4 l 5 l 6 l 36 l 37 l 38 l 



1=0 



2tt 



1 ho 
2 



(1) 



x 



{ln[l-r(jte,^)e~ 2 ^ 



hi 



1 



-2aqi 



where the reflection coefficients for two independent polarizations of electromagnetic 
field are defined as 



£iqi 



(2) 



£iqi + h ' h + qt 

Here k± is the magnitude of the wave vector in the plane of plates, = 2nkBTl/h 
are the Matsubara frequencies, ks is the Boltzmann constant, ei — e(i£z), and 



qi 



(3) 
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The problems in the application of the Lifshitz theory to real materials discussed 
in the Introduction are closely connected with the values of the reflection coefficients 
at zero Matsubara frequency. For later use we discuss it for the various cases. 

• For ideal metald^ it holds 

r||(0 > fcj.)=r±(0,*:x) = l. (4) 

• For real metals described by the dielectric function of the Drude model, 

wher e w p is the plasma frequency and v(T) is the relaxation parameter, it 



holds™ 



r||(0,*:j.) = l, rj.(0,*:j.) = 0. (6) 



Eq. ([6]) results in the discontinuity between the cases of ideal and real metals 
and leads to the violation of the Nernst heat theorem for metallic plates 
having perfect crystal latticeslS^EjUSHSEl 

• For real metals described by the dielectric function of the plasma model, 

em = l + % (7) 

from Eq. © it 

follows™] 

Y c?k\ + cj| - cfc ± 

r,| (0, fcj.) = 1, r ± (0,k±) = ^ r . (8) 

Jc?k\ + + ck ± 

Here, in the limit of ideal metals (u> p — > oo) the continuity is preserved 
because r±(0,k±) in Eq. © goes to unity. The free energy {TJ calculated 
with the permittivity ((5J) is also consistent with thermodynamics. 

• For dielectrics and semiconductors the dielectric permittivities at the 
imaginary M atsubar a frequencies are given by the Ninham-Parsegian 
representationp21Z0l 

j 3 

where the parameters Cj are the absorption strengths satisfying the condi- 
tion 

Y^Cj = £o -l (10) 

3 

and ujj are the characteristic absorption frequencies. Here, the static dielec- 
tric permittivity e = e(0) is supposed to be finite. Although Eq. (fT0|) is an 
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approximate one, it gives a very accurate description for many materials!^ 
By the substitution of Eq. ([9|) in Eq. ([2]) one arrives at 

r ll (0 ) k x ) = r o = ^^, rj.(0,A x ) = 0. (11) 

Note that the vanishing of the transverse reflection coefficient for dielectrics at zero 
frequency in Eq. (jll[) has another meaning than for the Drude metals in Eq. ([6]). 
For Drude metal the parallel reflection coefficient is equal to the physical value for 
real photons at normal incidence, i.e., to unity, and the transverse one vanishes 
instead of taking unity, its physical value. This results in the violation of the Nernst 
heat theorem for perfect crystal lattices. In the case of dielectrics both reflection 
coefficients at zero frequency in Eq. Ijlip depart from the physical value for real 
photons which is equal to ( v /eo — l)/(y / £o + 1). In this case, however, one of them is 
larger and the other one is smaller than the physical value. As we will see below, this 
leads to the preservation of Nernst 's heat theorem confirming that Eq. ((9]), despite 
being approximate, describes the material properties of dielectric and semiconductor 
plates in a thermodynamic consistent way. 

Now we derive the analytic representation for the Casimir free energy in Eq. (fT]) 
at low temperatures. For convenience in calculations, we introduce the dimcnsionlcss 
variables 

c & = ^ = Ti y = 2aqi , (12) 

where £ c = c/(2a) is the characteristic frequency, r = AirkBaT/ (he), and qi was 
defined in Eq. ([3]). Then the Lifshitz formula ((T|) takes the form 



3C 



where 

f(Cy) = f\\(C,y) + f±(Cv), (14) 

/||a(C,y)=yln[l-rf| a (C,2;) e - y ] , (15) 
and reflection coefficients in terms of variables (fl"2"|) . being given by 



eiy + ^y 2 + (f{ei~i) Vy 2 + Cf(ei-i) + y 

To separate the temperature independent contribution and thermal correction 
in Eq. (fT3|) we apply the Abel-Plana formulapEl 



Ul -%) m .r m « +i r«wk ! z&, (17) 



2 

1=0 v 

where -F(^) is an analytic function in the right half-plane. Here, taking it as 

/>OC 

F(x)= / dyf(x,y) (18) 
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and using Eq. (|17[) . we can identically rearrange Eq. (|13[) to the form 

T{a, T) = E(a) + AT (a, T) , (19) 

where E{a) is the energy of the van der Waals or Casimir interaction at zero tem- 
perature, 

and AT {a, T) is the thermal correction to this energy, 

. , . ihcr f 00 , F(irt) - Fi-irt) , . 

A ^ = ^j * K (21) 

Note that, in fact, Eq. (j2Tj) describes the dependence of the free energy on the tem- 
perature arising from the dependence on temperature of the Matsubara frequencies. 
Thus, AT(a, T) in (f2"Tj) coincides with the thermal correction to the energy, defined 
as T(a,T) — JF(a,0), only for plate materials with temperature independent prop- 
erties. 

The asymptotic exp ressions for the energy E{a) at both short and large sep- 
arations are well known jSEZEHl Below we find the asymptotic expressions for the 
thermal correction ([2~l|) under the conditions t«1 and r> 1. Taking into account 
the definition of r in Eq. (fT2|) . the asymptotic expressions at r < 1 are applicable 
both at small and large separations if the temperature is sufficiently low. 

We begin with condition r <C 1. Let us substitute Eq. ([9|) in Eqs. (fl4|) - (fTBj) . 
expand the function f(x,y) in powers of x = rt, and than integrate the obtained 
expansion with respect to y from x to infinity in order to find F(x) in Eq. (|18j) and 
F{ix) - F(-ix) in Eq. (j2T|) . 

It is easy to check that /x(C> V) does not contribute to the leading, second, order 
in the expansion of F(ix) — F(—ix) in powers of x. Thus, we can restrict ourselves 
by the consideration of the expansion 

/,, (x, y)=y]n(l- rge"") + , - ? , + , ^ ? + 0(z 3 ), (22) 

where ro was defined in Eq. (jlip . Note that for simplicity we consider here only one 
oscillator in Eq. © and put uij — u>\. The case of several oscillator modes can be 
considered in an analogous way. 

As a next step, we integrate Eq. (|22|) term by term according to Eq. (| 18[) . expand 
the partial results in powers of x and sum up the obtained series. Thereby we obtain 
the following expressions: 



Z 1 (x)= / ydyln(l-r 2 e-y)=-J2 



(1 + nx)e 



,,2» 



^(r^-^lna-r^ + O^ 3 ), (23) 
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Z 2 (x) = 



Z 3 (x) 



so + 1 

4rg^ac 



dy 



V (e y - rg) 
ydy 



2eqx 

So + 



T ^r 2 »Ei(-nx) 



(24) 



(1 + nx)e 



{s + J x ev - r 2 ( £o + ^ 
4£ c 2 



a^LiaW) 



2 1-rg 



+ 0(x 5 ) 



(25) 



(e + l)u; 2 

where Li„(^:) is the polylogarithm function and Ei(z) is the exponential integral 
function. 

From these equations it follows 

Zi(ix) - Zi{-ix) = 0(x 3 ), 



Zi{ix) — Z2{—ix) = 2i7r ■ 



so 



1 1 



Z 3 {ix) - Z 3 (-ix) = 0(x 5 ) , 
■ x 2 + 0(x 3 ) , 



(26) 
(27) 



and, thus, Z\ and Z 3 do not contribute to the leading order in the expansion of 
F(ix) — F(—ix). The latter is determined by Zi only. As a result, we arrive at 



F{ix) - F(-ix) = 



( £ o-l) 2 



iacc 3 + 0(a; 4 ), 



(28) 



2(e + l)' 

where ro was substituted from Eq. and a was introduced for the still unknown 
real coefficient of the next to leading order resulting from Z± and Z2 as well as, 
possibly, from f±(C,y)- At this stage it is difficult to determine the value of this 
coefficient because all powers in the expansion of f(x,y) contribute to it. Remark- 
ably, the two leading orders depend only on the static dielectric permittivity So and 
are not influenced by the dependence of the dielectric permittivity on the frequency 
contained in Z 3 . 

Substituting Eq. (28|) in Eq. (JH) and using Eq. fl9j). we obtain 



T(a,T) = E(a) 



he 



327r 2 a 3 



C(3) ( £ Q-1) 2 3 

8tt 2 £ + 1 



C4T + 0(r 5 



(29) 



where C4 = a/240 and £(z) is the Riemann zeta function. 

So far we have considered the free energy. The thermal pressure is obtained as 

he 



P(a,T) = -^g^ = Po(a) - [C 4 r 4 + 0(r 5 )] 



(30) 



9a ~ uv " y 327r 2 a 4 

where -fo( a ) — —dE(a)/da is the Casimir pressure at zero temperature. 

In order to determine the value of the coefficient C4 of the leading term, we 
express the pressure directly through the Lifshitz formula 



P(a,T) 



her 

~327r 2 a z 



;=n 



(0,y) 



ri(0,y) 



e«-r 2 (0,y) e» - r 2 (0 , y) _ 



(31) 



Again, applying the Abel-Plana formula (|17p . we represent the pressure as follows, 
P(a,T) =Po(a) + AP(a,T), (32) 
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where the thermal correction to Pq(ci), the pressure at zero temperature, is 

ihcr 



AP(a,T) = 



327r 2 a 4 

and the function <I>(x) is given by 



00 $(iri) - $(-iTt) 



- 1 



$(x) = $|| (x) + $j_(a 



y 2 rf l± (x,y) 

e» - 



0,2/) 



(33) 



(34) 



First, we determine the leading term of the expansion of $j_(x) in powers of 
x. For this purpose, let us introduce the new variable v — y/x and note that the 
reflection coefficient r±(x,v) depends on x only through the frequency dependence 
of e given by Eq. ©. Thus, we can rewrite and expand Eq. (|3"4"|) as follows: 

v 2 r\[x, v) 



$_l(x) = x i J dv ■ 
where, according to Eq. H 



r\{x,v) 



dv ■ 



1 



l(v) 



0(x q 



rj_(v) = rx(0,v) = 



V^ 2 



\/v 2 + £o - 1 + V 

Integration in Eq. (j35|) with account of Eq. (|36|) results in 



e (3- e ) 



from which it follows: 



$±(ix) - $±(—ix) 



2^ 



1 



6 



+ 0(x 4 ), 



(35) 



(36) 



(37) 



(38) 



The expansion of $11 (x) from Eq. p4p in powers of x is somewhat more cumber- 
some. It can be performed in the following way. As is seen from the second equality 
in Eq. (|26p . the dependence of the dielectric permittivity on frequency contributes 
to F(ix) — F(—ix) starting from only the 5th power in x. Bearing in mind the con- 
nection between free energy and pressure, we can conclude that the dependence on 
the frequency contributes to ^ii(ix) — ix) starting from the 4th order. We are 
looking for the lowest (third) order expansion term of <&ii (ix) — <J>ii (— ix). Because of 
this, it is permissible to disregard the frequency dependence of e and describe the 
dielectric by its static dielectric permittivity. 

To begin with, we identically rearrange $11 (x) in Eq. (f34)) by subtracting and 
adding the two first expansion terms of the function under the integral in powers 
of x, 



$|| (x) 



dy 



y 



r jjO,y) 



2e n 



y 



£o + l(l 



r 2 e-v) 



+ 



y 2 dy 



2er 



ey 



£0 + 1 



dy 



(1 - r 2 e-vy 



(39) 
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and consider these three integrals separately. The first integral in terms of the new 
variable v — y/x reads 



Qi(x) = x d 



dv 



2e 
£o + 1 (1 



r 2 -vx 

r e 



(40) 



s v - ^v 2 +s a -l 



(41) 



where, in accordance with Eq. (I1G|) . 

= r \\(°, v ) -> r 

EqV + v?r + so — 1 

Expanding Qil^O m powers of x and explicitly calculating the remaining inte- 
grals for the lowest, third, power of x results in 



Qi(x) = x' 



dv 



1-rKv) 



2e n 



1 - rg e + 1 (l _ r 2 ) 2 



~6~ 



l-^(l + 3 £o -2 £ 2 ) 



2r 2 



12 £ 



+ 0(.t 4 ) (42) 



+ 0(x 4 ). 



l-»o £ o + l(l-r 2 ) 2 . 

The second and third integrals on the right-hand side of Eq. (|39[) are simply deter- 
mined with the following result: 



Q 2 (x) 



3Li 3 (r 2 ) - 

2e a; 2 
eo + 1 

2e 



2 „ (2 + 2nx + n 2 a; 2 )e 



2_2\„ — in 



3 1 



^ + 0(x 4 ). 



2e a; 2 



{l-rle-yf 



eo + 1 



1 



1 



eo + 1 1 - rle~ x 
+ 0(x 4 ). 



'0 \ x _ 'OJ 

Substituting Eqs. [|l2]). (|43jl and (44]) into $y(a;) = Qi(a;) + Q 2 (x) 
arrive at 



$H (ix) - $|| (—ix) = -i 



1 



2e 2 - 3e - 1) 



+ 0{x 4 ). 



Then, by summing Eqs. 



and (|45p. the result is obtained 



$(ix) - $(-ws) = -i ( 2 + £o /2 - £o /2 - 2 V^) V + °( x 



(43) 

(44) 
? 3 (i), we 

(45) 
(46) 



Now we substitute Eq. (|46p in Eq. (|33|) and perform integration. Finally, from 
Eq. (|32p the desired expression for the Casimir pressure is derived 



P(a,T)=P (a) 



he 



32ir 2 a 4 



(V^-l) + r4 , Q(r5) 



720 



(47) 



By comparison with Eq. (|30|) the explicit form of the coefficient C4 is found as 

Ci= (v^-l)(£o + £o /2 - 2)/720 (48) 
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and, thus, both two first perturbation orders in the expansion for the free energy 
(|29| are determined. 

Equations ([29|) , (|47| and (j48| solve the fundamental problem of the thermody- 
namic consistency of the Lifshitz theory in the case of two dielectric plates. From 
Eqs. (|29|) and (j48|) the entropy of the van der Waals and Casimir interaction between 
plates takes the form 



S(a,T) = 



dF(a,T) 3fc s C(3)(eo - I) 2 



dT 647r 3 a 2 (e + 1) 

2ir 2 (eg + 1) {el 12 + 2e + 2^ + 2) 
135C(3)(V^+1) 2 



(49) 



r + 0(r 2 ) 



As is seen from Eq. (l4"9")l , in the limit r — * (T — * 0) the lower order contributions 
to the entropy are of the second and the third powers in the small parameter r. 
Thus, the entropy vanishes when the temperature goes to zero as it must be in 
accordance with the third law of thermodynamics (t he Xernst heat theorem) . 

A similar behavior was obtained for ideal metal J 4 " ! 72 ' 7 -^ and for real metals 
described by the plasma model P^^ For example, in the case of plates made of 
ideal metal the entropy at low temperatures is given bjj4QlZll73] 



S(a,T) 



3fcaC(3) 
327r 3 a 2 



1 



2^ 2 



135C(3) 



0(r 2 ) 



(50) 



Note, however, that the expansion coefficients in Eq. ([50]) cannot be obtained as 
a straightforward limit |eo| ~ > 00 in Eq- (|4"5|) and the above equations for the 
free energy and pressure. The mathematical reason is that it is impermissible to 
interchange the limiting transitions r — > and |eo| 00 i n the power expansions 
of functions depending on £o as a parameter. 

Remarkably, the low-temperature behavior of the free energy, pressure and en- 
tropy of nonpolar dielectrics in Eqs. ([29]) . (|47|) and (|49)l is universal, i.e., is deter- 
mined only by the static dielectric permittivity. The absorption bands included in 
Eq. ([9]) do not influence the low-temperature behavior. A more simple derivation 
of the results ([29]) . (|47|) (|49|) for dielectrics with constant e is contained in Ref. [74l 
As was demonstrated above, all these results remain unchanged if the dependence 
of dielectric permittivity on frequency is taken into account. 

In Ref. [60] more general results were obtained related to two dissimilar dielectric 
plates with dielectric permittivities e^(ia) and e^ 2 ^(w). For brevity here we present 
only the final expressions for the low-temperature behavior of the Casimir free 
energy, pressure and entropy between dissimilar plates. They are as follows!^ 



F{a,T) = E{a) 



he 



327r 2 a 3 



03) 
8tt 2 



.(2) 



+ 2e. 



aU2) (jx) 



c 



- D(4 2) - 1 



(4 X) 



D(4 2) 



.(2h 



C 4 r + 0(r 2 ) 
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P(a,T) = P (a) 



S(a,T) = 



he 



32ttV 



[C 4 r 4 + 0(r 5 )] 



(51) 



2na 2 



3C(3) 

32. 2 ( £ « + l)(4 2 ) + l) ( 4 1} + 4 2) ) 

Here Eq 1 ' 2 ^ = e( 1,2 )(0) and the coefficient C4 is given bjJ^ 



-C 4 r + 0(r 2 ) 




(52) 



It is easily seen that in the limit Eq 



(i) 



.(2) 



eo equations (|51j) . ([52| coincide with 



equations (|29p . ((171) - P5)) having obtained above. Note that in the application region 
of low-temperature asymptotic expressions the entropy of the Casimir interaction 
between dielectric plates is nonnegative. 

The obtained analytic behavior of the free energy, pressure and entropy at low 
temperatures can be compared with the results of numerical computations using 
the Lifshitz formula. Dielectric properties of the plates can be described by the 
static dielectric permittivity or more precisely using the optical tabulated data for 
the complex index of refraction. As an example, in Fig. 1 we present the thermal 
corrections to the Casimir energy (a) and pressure (b) at a separation a = 400 nm 
as functions of temperature in the configuration of two dissimilar plates made of 
high- resistivity Si and Si02- The dielectric permittivities of both materials along 
the imaginary frequency axis were computed in Ref. [75] using the optical data of 
Ref. [75] The precise thermal corrections computed by taking into account these 
permittivities are shown by the solid lines and corrections computed by our analyt- 
ical asymptotic expressions are shown by the long-dashed lines. Short-dashed lines 
indicate the results computed by the Lifshitz formula with constant dielectric per- 



mittivities of Si and Si02 equal to Eg = 11.67 and e^' = 3.84, respectively. As is 
seen in Fig. la,b, at T < 60 K the results obtained using the analytical asymptotic 



(2) 
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-A^(pJ) 




25 50 75 100 125 150 175 200 




25 50 75 100 125 150 175 200 

T(K) 

(b) 



Fig. 1. Magnitudes of the thermal corrections to the energy (a) and pressure (b) in configuration 
of two plates, one made of Si and another one of Si02, at a separation a = 400 nm as a function of 
temperature calculated by the use of different approaches: by the Lifshitz formula and tabulated 
optical data (solid lines), by the Lifshitz formula and static dielectric permittivities (short-dashed 
lines), by the asymptotic expressions in Eqs. I|51| l and (l-">_'l (long-dashed lines). 



expressions practically coincide with the solid lines computed using the tabulated 
optical data for the materials of the plates. 

Now we return to the case of two similar dielectric plates and consider the 
asymptotic expressions unde r the con dition r ^> 1, i.e., at high temperatures (large 
separations). It is well that in this case the approximation of static 
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dielectric permittivity works good and the main contribution is given by the zero- 
frequency term of the Lifshitz formula ()13|) 



ftrr f°° 

^ (a ' T)= 6^W V*vHl-&^) (53) 



(the other terms being exponentially small). Performing the integration in Eq. (|53j) 
we obtain 

T ^ T) = -^ U3 ^- (54) 
In a similar manner for the Casimir pressure and entropy at t 3> 1 it follows 

P(a ' T) = "8§ Li3 ™ ' S{a ' T) = 1^ Li3 <® ■ (55) 
Equations ([54]) and ([55]) are simply generalized^ for the case of two dissimilar di- 
electric plates by performing the replacement rg — ► r^Tg , where ?"q are defined 
by Eq. (fTT|) with the static dielectric permittivities of dissimilar plates e^' 2 \ 

3. IS THE DC CONDUCTIVITY RELATED TO THE CASIMIR 
INTERACTION BETWEEN DIELECTRICS? 

As was discussed in the previous section, the zero- frequency term in formula (|13p . 
i.e., the contribution with 1 = 0, is of prime importance and determines many of 
the basic properties of the Casimir interaction. In the above consideration we have 
described dielectric materials by Eq. ((9]) with finite static dielectric permittivity £o- 
This resulted in Eq. where one reflection coefficient at zero frequency is larger 
and the other one is smaller than the physical value for real photons at normal 
incidence. However, in the Lifshitz theory, the departure of both coefficients from 
their physical values is coordinated in such a way that the Nernst heat theorem 
remains valid. 

It is common knowledge that at nonzero temperatures dielectric materials pos- 
sess a negligibly small but not equal to zero dc conductivity. From physical intuition 
it is reasonable to expect that the influence of this conductivity on the van der Waals 
and Casimir forces should be also negligible. In Ref. [57] it was shown that, on the 
contrary, the inclusion of small dielectric dc conductivity in the model of dielectric 
response leads to a large effect in dispersion forces. This raises the question if dc 
conductivity is related to dispersion forces or if, on the contrary, the zero-frequency 
contribution should be understood not literally but as an analytic continuation from 
the region of high frequencies determining the physical phenomenon of dispersion 
forces. 

To illustrate this problem in more details, we consider the asymptotic behavior 
of the free energy and entropy at low temperature with included dc conductivity. 
What t his mea ns is that, instead of the dielectric permittivity e(i£i) given in Eq. ([9]), 
one use d 57 * 58 ! 

mi)=^m + ^=^l(iO + ^p-, (56) 

6 1 



February 5, 2008 2:30 WSPC/INSTRUCTION FILE tcxtDr 



16 B. Geyer, G. L. Klimchitskaya and V. M. Mostepanenko 



where oo is the dc conductivity of the plate material and (3{T) = 2hao/(k B T). The 
conductivity of dielectrics depends on temperature as <tq ~ exp(— b/T) where b is 
determined by the energy gap A which differs for different materials!^! The small- 
ness of the dc conductivity of dielectrics can be illustrated^ by the example of 
Si0 2 where at T = 300 K it holds (3 ~ 10" 12 . Thus, the role of the dc conductivity 
is really negligible for all I > 1. In addition, (3{T) quickly decreases with decrease 
of T and, as a consequence, remains negligible at any T . In spite of this, the substi- 
tution of Eq. (|56|) into the reflection coefficients (fl6|) leads to different result than 
in Eq. ([IT]): 



f||(0,y) = l, f||(Q,y) = 0. (57) 

Equation (|57|) is in some analogy to Eq. ([6]), obtained for metals described by the 
Drude model, which leads to the violation of the Nernst heat theorem in the case 
of perfect crystal lattices. 

Now we substitute the dielectric permittivity £/ = in Eqs. (fT3| - (fT6|) 

instead of Ei and find the Casimir free energy J- (a, T) with included dc conductivity. 
For convenience, we separate the zero-frequency term, subtract and add the usual 
zero-frequency contribution for dielectric without dc conductivity. The result is 
k B T 



T{a,T) 



+ 



+ 



167TO 2 

k B T 
167ra 2 

k B T 



ydy 



In 1 - 



-In 1- 



ydy In (l 



(58) 



°° f°° r r i r i -i 

J2 / ydy {in l-rf(Q,y)e-y +ln 1 - r~i(0, y)e- y } 



1=1 "'w 



where tq was defined in Eq. (jll|) . Let us expand the last, third, integral on the right- 
hand side of Eq. (|58[) in powers of the small parameters (3/1. Then, we combine the 
zero-order contribution in this expansion with the second integral on the right- 
hand side of Eq. ([58| and obtain the Casimir free energy f(a, T) calculated with 
the dielectric permittivities ei . The first integral on the right-hand side of Eq. (1ST 
is calculated explicitly. Then, Eq. (|58|) can be rewritten as 

k B T 



F(a,T) = T(a, T) - 



[C(3)-Lia {rl)\ +R(a,T). 



(59) 



lQna 2 

Here, R(a,T) is of order 0((3/l). It represents the first and higher-order contribu- 
tions in the expansion of the third integral on the right-hand side of Eq. (|58|) hi 
powers of (3/1. Restricting its explicit form to the first order contribution we get 



R(a,T) = Ri(a,T) 
k B T 



Ri(a,T) 



OO 

\ ^ 



0[{(3/lf 



ft 



Ana 2 ^ I 
i=i 

r \\(Chy) 



dy y 2 e~ 



(2-e l )C 2 -2y 2 



Cl ^y 2 + Q{e l -l) 



[Vy 2 + Cf(ei-i) 



l-r 2 (Ci,y)e-y y y 2 + C 2 {£l _ 1)+ , 



(60) 



(61) 
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Calculating the entropy by the first equality in Eq. (|49|) . we arrive at 

S(a,T) = S(a,T) + [ C (3) - Li 3 (r 2 )] - (62) 

where S(a, T) is the entropy for the plates with the dielectric permittivity e; given 
by Eq. (|49]). 

Let us now show that the quantity R(a, T) exponentially goes to zero with the 
decrease of T. First we consider only the integral in Eq. (j6Tj) , expand the integrated 
function in powers of r (we recall that Q = tI), restrict ourselves to the main 
contribution, resulting for r = 0, and rearrange it appropriately: 

, ye- y r 2 f°° r 2 a e-y 

ay (e + ir l-rle-v (e 2 -l)J Cl ^l-r^ 

V I n —i JQ \ I n —x 

Substituting this in Eq. (|61|) . we find 

v u ; n=l \l=l 1=1 



e -nrl 



in ,,- ( T ; - 1) „ 



2 _ 1-V r,2 



ln(l 



fcfl Lio (Vn ) , n* 



(64) 



Here, the last line is obtained by using the equality 

-ln(l-e-" T ) + UT = -lnT + l-lnn + 0(r 2 ), (65) 
substituting only its leading term and observing the definition of the integral log- 

oo 

arithm, Li 2 (z) = (1/2) J2 z n /n 2 . Taking into account that (3 ~ (1/T) exp(-b/T), 

n=l 

we get the conclusion that the temperature dependence of R±(a,T) is given by 

Ri(a,T) ~ e- b/T \nT. (66) 

Thus, both Ri(a,T) and its derivative with respect to T in Eqs. ([55]) and (jl?2")) go 
to zero. The terms of the second and higher powers in (3 in Eq. (|60[) go to zero even 
faster than Ri when T — > 0. 

Finally, in the limit T — > from Eq. (|62|) it follows 



^(«,0) = I ^[C(3)-Li 3 (r 2 )] >0. (67) 

The right-hand side of this equation depends on the parameter of the system under 
consideration (the separation distance a) and implies a violation of the Nernst heat 
theorem. An analogous result was obtainecP^in the case of two dissimilar dielectrics. 
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The violation of the Nernst heat theorem in the Casimir interaction for dielectrics 
originates from the inclusion of the dc conductivity in the model of dielectric re- 
sponse. This violation is, however, of different nature than the one discussed above 
in the case of Drude metals. In the case of dielectrics the entropy at zero tempera- 
ture is positive but in the case of Drude metals it is negative. In the case of metals 
the violation is caused by the vanishing contribution from the transverse electric 
mode at zero frequency whereas the other reflection coefficient takes the physical 
value 1 [see Eq. ([£])]. For dielectrics the situation is quite opposite. In this case the 
transverse reflection coefficient at zero frequency is always equal to zero [compare 
Eqs. (fTT|) and ([57| in the absence and in the presence of the contribution from 
dc conductivity]. Here the violation occurs due to the unity value of the parallel 
reflection coefficient at zero frequency in Eqs. (|57|) which departs from the value 
? 'o = (^o ~ l)/( £ o + 1) coordinated with the zero value of the transverse coefficient 
in Eq. (fTTjl. 

One can conclude that the dc conductivity of a dielectric is not related to the na- 
ture of the van der Waals and Casimir forces and must not be included in the model 
of dielectric response. Ignoring this rule results in a violation of thermodynamics. 
Physically it is amply clear that there is no fluctuating field of zero frequency and 
that for such high-frequency phenomena as the van der Waals and Casimir forces the 
low-frequency behavior should be obtained by analytic continuation from the region 
of high frequencies. This permits to conclude that the correct procedure consists in 
the substitution of the finite static dielectric permittivities into the zer o-frequen cy 
term of the Lifshitz formula, as Lifshitz and his collaborators really didJSSEZEEl 

4. THERMAL CASIMIR FORCE BETWEEN DIELECTRIC 
AND METAL PLATES 

The Casimir interaction between metal and dielectric plates suggests the interest- 
ing opportunity to verify the thermodynamic consistency of the Lifshitz theory 
with different models of the dielectric response. This configuration was first inves- 
tigated in Ref. [61] where it was proved that the Casimir entropy is in accordance 
with the demands of the Nernst heat theorem if the static permittivity of the di- 
electric plate is finite. In Ref. I61[ however, only the first leading terms in the low- 
temperature asymptotic expressions for the free energy and entropy were obtained 
and the Casimir pressure was derived only in the dilute approximation. Here we 
derive the more precise low-temperature behavior for the Casimir free energy, pres- 
sure and entropy in the configuration of one plate made of ideal metal and another 
plate made of dielectric with any finite static dielectric permittivity. 

For the configuration of metal and dielectric plates the Lifshitz formula takes 
the forir|2QI analogical to Eq. (fl~3"|) 




(68) 
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x {in [l-rf (Ci,ff)rf (Ci,tf)e-*] + In [l - rf (£,, y) rf (0, y) e~»] } . 

Here the reflection coefficients ^^j^ for metal and dielectric, respectively, are given 
by Eqs. IT6|) where e; should be changed for e\ l ' D — E M ' D (i£i). 

For an ideal metal r^ ± (Ci,y) = 1 and Eq. (|68|) takes the more simple form 

*^-sH?i:(i-|)/> <«» 

x {ln[l-r,|(Ci,tf)e-*] + In [l - rxfe,^] } 

(here and below we omit the index D near the reflection coefficient and permit- 
tivity of a dielectric plate). We admit that the dielectric permittivity calculated at 
Matsubara frequencies £; = £o, i-e-j is equal to its static value and find the asymp- 
totic behavior of Eq. (|69|) at small r. [In analogy with Sec. 2 it is possible to prove 
that the deviations of e(i£;) from Eq at high frequencies do not influence the low- 
temperature behavior of the Casimir free energy, pressure and entropy. It can be 
shown also that the results of this section are valid not only for ideal metal plate 
but for plate made of real metal as well.] The free energy (|69|) can be represented 
by Eqs. fTSj) - (jUJ), with the function /(C,y) replaced by 

/(£, y) = y In [1 - r,| (C, y)e-y] + y In [l - rj_(C, y)e~ y ] 

= f\\(C,y) + f±(C,y)- (70) 

In the case of one dielectric and one metal plate both /|| and f± contribute to 
F(ix) — F(—ix). The expansion of f(x,y) in powers of x takes the form 



fo^ e -s/_^^ y>j- W2 + 0(z 3 ). (71) 

Now we integrate Eq. (|7T|) in accordance with Eq. §TE\i to find the function F(x). 
The integral of the first term on the right-hand side of Eq. (f7Tj) is evaluated using 
the new variable v = y — x: 

/>oo 

ydy\n(l-r er y ) = vdv ln(l - r^r v ) + 0{x 2 ), (72) 
Jo 

where the coefficient near the first-order contribution in x vanishes. As a result, this 
term could contribute to F{ix) — F(—ix) only starting from the third expansion 
order. The integrals of the second-order terms on the right-hand side of Eq. ([7T|) 
are simply calculated using the formulas 

dy = -Ei(-ar), / dy = -Ei(-nx). (73) 

y Jx y 

Finally, we obtain 

F(ix) - F(-ix) = m (£ ° ~ 1)2 N x 2 -i 1X 3 + 0(x 4 ), (74) 

4 (£ + l) 
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where the unknown third order expansion coefficient is designated as 7. 

Substituting Eq. ([74]) in Eq. (j21| and using Eq. (fl"9|) . we find the free energy in 
the system metal-dielectric in the form 



T(a,T) = E{a) 



he 



327r 2 a 3 



l07H £0 + 1 



(75) 



where if 4 = 7/240. 

The Casimir pressure in the configuration of metal and dielectric plates obtained 
from Eq. (|75p is equal to 

he 



P(a,T)=P (a) 



(76) 



327r 2 a 4 

The direct application of the Lifshitz formula gives the expression for the pressure 
analogical to Eq. ([BT]) . 



P(a,T) 



her 

327r 2 a 4 



E 



1=0 



2 



y 2 dy 



(ti,y) 



^ r±(Ci,y) 



(77) 



Using the Abel-Plana formula (fTT)) . Eq. (|77|) can be represented in the form of 
Eqs. AMI, (ESI) where 

y 2 r\\ y ±(x,y) 



dy- 



(78) 



eV - r iA x ,y) ' 

Again, we deal first with <fr±(x). By adding and subtracting the asymptotic 
behavior of the integrated function at small x, 



? r±iX : V \ =\(eo-lWe-v + 0(x>).. 
e v -r±(x,y) 4 



(79) 



and introducing the new variable v = y/x, the function $>±(x) can be identically 
rearranged and expanded in powers of x as follows: 



<Mx) = i(e -l> 2 e-* 



1 



-(e - l)ar(l - x) + x' 



dv v 2 r" (v)e~ 

- n=l 

v 2 r±(v) s -l 



- 2 ( £o -l)e— 



dv 



1 - rx(u) 



+ 0(.t 4 ). (80) 



The integral on the right-hand side of Eq. ([8H|) is converging and can be simply 
calculated with the result 



^x 2 - l -(e Q ^T a -l)x* + 0{x A ). 



(81) 



To deal with $11 (x) we add and subtract in Eq. ([75)1 the two first expansion 
terms of the integrated function in powers of x, 



$11 (*) 



y 2 dy 



e r e s x 2 



e y - r y i(e a + 1) (1 - r e~v) 
r \\(x,y) r 



(82) 



e r e y x 2 



ey~r\\{x,y) eV - r y 2 (e + 1) {1 - r Q e-yy 
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The asymptotic expansion of the first integral on the right-hand side of Eq. (j82[) is 
given by 



2Li 3 (r ) 



£o(eo - 1) 



2(eo- 

and of the second one by 



-L( £o -1)(3 £q -2)x 3 



By summing Eqs. (|83f and (f84|) we find 



0(x 4 



x 3 + 0(x 4 ). 



$11 0) 



2Li 3 (r ) 
1 



e Q {e - 1) 



2(eo + 1) 

j [(e - 1) + (eoV^ - !) - ^o(e a - 1)V^] x 3 + 0(x 4 ). 
Finally we add Eq. I 



(83) 



(84) 



(85) 



to Eq. I[8ip and arrive at 
2 



*3 (l"2e V^ + £o\/^)^ 3 + 0( a;4 )- 



(86) 



Substituting this in Eq. (|33|) and using Eq. (|32[) . the asymptotic expression for the 
Casimir pressure is obtained 



P{a,T) = P (a) 



he 



3^0 (1 - 2eoV^ + £ov^) r A + 0(r 5 ) 



(87) 



32?r 2 a 4 

Comparing Eqs. (|87|) and (|76|) . the explicit expression for the coefficient K4 reads 

^=3^(l-2e ov ^ + ^Ve5)- ( 88 ) 

Now we are in a position to find the asymptotic behavior of the entropy in the 
limit of low temperatures in the configuration of two parallel plates one of which 
is metallic and the other one dielectric. By calculating the negative derivative of 
Eq. I|75p with respect to temperature, one arrives at 



S(a,T) = 



3fc£C(3)(£o-l) 2 2 
1287r 3 a 2 (e + l) T 

8n 2 (e + 1)(1- 2e V^ + £ov / ^) 



(89) 



r + 0(r 2 ) 



135C(3)(e -l) 2 

This equation is in analogy to Eq. (j49|) obtained for the case of two dielectrics. 
As is seen from Eq. (|89|) , the entropy of the Casimir and van der Waals interactions 
between metal and dielectric plates vanishes when the temperature goes to zero, 
i.e., the Nernst heat theorem is satisfied. Note that the first term of order r 2 on the 
right-hand side of Eq. (f59"| was already obtained in Ref. [STJ It is notable also that at 
low temperatures the entropy goes to zero remaining positive. At the same time, as 
was shown in Ref. 1611 at larger temperatures entropy is nonmonotonous and may 
take negative values. This interesting property distinguishes the configuration of 
metal and dielectric plates from two dielectric plates. In the latter configuration the 
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negative entropy appears only for nonphysical die lect rics with anomalously large 
and frequency independent dielectric permittivities 

Now let us consider the case r ^> 1, i.e., high temperatures (large separations). 
In the same way as for two dielectric plates, here the main contribution to the free 
energy is given by the term with I = in Eq. (|69l) . 

Yict r°° 

Ha ' T)= 64^J ^ Ml -roe"*). (90) 
Performing the integration in Eq. (|90p , we obtain 

^ (a ' r)= -iS Li3(ro) - (9i) 

From this equation, for the Casimir pressure and entropy at r 3> 1 it follows 

P(a, T) = -g^Li 3 M , S(a, T) = ^U 3 (r ) . (92) 
The results (|9Tj) and (|92[) are analogous to (|54)l and ([55)) for two dielectric plates. 



5. THE PROBLEM ORIGINATING FROM THE ACCOUNT OF 
DIELECTRIC DC CONDUCTIVITY 

In the previous section it was supposed that the static dielectric permittivity £o of 
the dielectric plate is finite. Now we will deal with the configuration of metal and 
dielectric plates with included dc conductivity of the dielectric material. In doing 
so the permittivity of the dielectric plate is given by 

£(*&)= £0 + ^, (93) 

where all notations were introduced in and below Eq. (|56p . Thus, the reflection 
coefficients at zero frequency satisfy Eq. Q for a plate made of ideal metal and 
Eq. (|F7|) for a plate made of dielectric with included dc conductivity. Let us find 
the low-temperature behavior of the Casimir entropy and verify the consistency of 
the Lifshitz theory with thermodynamics in this nonstandard situation. 

For this purpose we substitute the dielectric permittivity (|93|) in Eq. (|69|) in- 
stead of £o and find the Casimir energy J 7 (a, T) with included dc conductivity of 
a dielectric plate. In the same way as in Sec. 3, it is convenient to separate the 
zero-frequency term of ~F(a, T) and subtract and add the usual zero-frequency con- 
tribution for metal-dielectric plates computed with the dielectric permittivity eq, 

j rp i>00 

+ £r^ rVvM 1 -roe'*) (94) 

h T /"CO 

+ ^E/ Vdy{hi[l-r\\(b,v)e-y]+]n[l-r x (Ci,v)e- v ]}. 
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Here the reflection coefficients fiu are calculated with the permittivity (|93"1) . We 
expand the third integral on the right-hand side of Eq. (|94|) in powers of the small 
parameter (3/1. The zero-order contribution in this expansion together with the 
second integral of Eq. (|94|) form the Casimir free energy T(a,T) calculated with 
dielectric permittivity Eq. The first integral on the right-hand side of Eq. (|94ll is 
calculated explicitly. As a result, Eq. (|94[) is rearranged to 

F(a, T) = F(a, T) - [C(3) - Li 3 (r )] + Q(a, T), (95) 

where Q(a,T) contains the first and higher-order contributions in the expansion of 
the third integral on the right-hand side of Eq. (|94p in powers of (3/1. The explicit 
form of the main first-order term in Q(a, T) is the following: 

n< T , k B T^(3 [°° dyy 2 e~y \ {2-e )Q-2y 



C. vV + Cf(eo -1)1 y v * + Cf(e —) + soy] 2 

x r—. .. ^ ~ 9 7^ \- (96) 

l-r|,(Ci,»)e-» y y 2 + C f {£o _ 1)+y Y l-r x (Ci,y)e-"J 

In the same way as in Sec. 3, we expand the integrated function in Eq. (|96p in 
powers of r (bearing in mind that Q = tI) and preserve only the main contribution 
at r = 0: 



l ( t\ k B T ^ (3 f°° dyyr e y 
hM = -^^li (el -l)(l-r e-) (9?) 



k B Tf3 



\ -> £o_ 



00 ~ — nrl x 

-i i=i 



E 



W( £ 2 - 1) ^ 

Dealing with this expression in the same way as with Eq. (|64p , we arrive at 

Qi(a,T) ~ e~ b/T lnT. (98) 

The Casimir entropy in the configuration of metal and dielectric plates with 
included dc conductivity of the dielectric plate is obtained as minus derivative of 
Eq. (|9"5|) with respect to temperature, 

kB fsfi\ t t m <9Q(a,T) 
W [C(3) - Ll3M] -^T- 
Using Eq. (|98|) . the calculation of the limiting value at T — > is straightforward: 

S(a,0) = -^[C(3)-Li3(r )]>0. (100) 

From this equation it follows that the inclusion of the dc conductivity of dielectric 
plate in the configuration metal-dielectric results in a violation of the Nernst heat 
theorem. In the above this result was obtained for a metallic plate made of ideal 
metal. It can be shown that it remains valid for a metallic plate made of real metal 
with finite conductivity. 



S(a, T) = S(a, T) + ^ [C(3) - Li 3 (r„)] - (99) 
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Thus, both configurations (two dielectric plates or one metal plate and one 
dielectric) lead to the same conclusion that when the dc conductivity is included 
in the model of dielectric response of the dielectric plate, the Lifshitz theory loses 
its consistency with thermodynamics. This confirms the conclusion made in Sec. 3 
that the actual properties of dielectric materials at very low frequencies are in fact 
not related to the van der Waals and Casimir forces. 

6. QUALITATIVE DISCUSSION OF CASIMIR INTERACTION 
BETWEEN METAL AND SEMICONDUCTOR 

As was mentioned in the Introduction, semiconductors possess a wide variety of 
electric and optical properties ranging from metallic to dielectric. This opens the 
possibility to modulate the van der Waals and Casimir forces by changing the charge 
carrier density. Bearing in mind the discussed above problems on the consistency 
of the Lifshitz theory with thermodynamics, semiconductors can provide us with 
a test for the validity of different approaches. Thus, if for good dielectric the dc 
conductivity does not play any real role in the van der Waals and Casimir forces, 
the question arises on how much it should be increased in order to become a relevant 
factor in the description of dispersion forces. 

In Ref. [8T]the Casimir force acting between two Si plates was calculated using the 
simple analytic expression for Si dielectric permittivity as a function of frequency. 
The complete tabulated optical data of Si were used in Ref. [75] to calculate the van 
der Waals interaction of different atoms with a Si wall. The first attempt to measure 
the van der Waals force between a glass lens and a Si plate and to modify it by light 
due to the change of carrier density was undertaken in Ref. 82) However, glass is a 
dielectric and therefore the electric forces due to localized point charges could not 
be controlled. This might explain that no force change occured on illumination at 
small separations where the effect should be most pronounced. 

The first measurements of the Casimir force between a gold coated sphere and 
a single crystal Si plate were performed in Refs. [62l [63] by means of the atomic 
force microscope. The experiments used a p-type B doped Si plate of resistivity 
p = 0.0035 Q cm. The chosen resistivity of the plate is in some sense intermediate 
between the resistivity of metals (which is usually two or three orders of magnitude 
lower) and the resistivity of dielectrics (it can be by about a factor of 10 5 larger; for 
instance, high-resistivity "dielectric" Si has the resistivity po = 1000 SI cm). Thus, 
the used Si plate had a relatively large absorption typical for semiconductors but it 
was also enough conductive to avoid the accumulation of charges!^ 

In Fig. 2, taken from Ref. [531 the differences of the theoretical and mean ex- 
perimental Casimir forces acting between Au sphere and Si plate are presented as 
functions of separation. In Fig. 2a the theoretical force i^ thcor is computed using 
the Lifshitz formula and the dielectric permittivity of a Si plate with the relatively 
low resistivity p used in experiment. This dielectric permittivity goes to infinity as 
with decreasing frequency (see the solid line in Fig. 3). In Fig. 2b the theo- 
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Fig. 2. Differences of the theoretical and mean experimental Casimir forces versus separation. 
Theoretical forces are computed (a) for the Si plate used in experiment and (b) for dielectric Si. 
Solid and dashed lines indicate 95 and 70% confidence intervals, respectively. 

retical Casimir force F theor is computed using the dielectric permittivity of the Si 
plate made of "dielectric" Si with high resistivity p . The dielectric permittivity of 
high-resistivity Si is shown by the dashed line in Fig. 3. It is characterized by the 
finite static value e Sl (0) = 11.67. The solid and dashed lines in Fig. 2a, b indicate 
the 95% and 70% confidence intervals, respectively. As is seen from Fig. 2, the the- 
oretical approach using the dielectric permittivity of high-resistivity "dielectric" Si 
is excluded by experiment within the separation range from 60 to 110 nm at 70% 
confidence. At the same time, the theory using the dielectric permittivity of Si with 
a low resistivity p is consistent with experiment. 
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Fig. 3. Dielectric permittivity of Si plate used in experiment along the imaginary frequency axis 
(solid line). Dashed line shows the dielectric permittivity of dielectric Si. 



The above results suggest an approach on how to correctly determine the pos- 
sible role of the low-frequency conductivity properties in dispersion forces. As is 
seen from Fig. 3, the dielectric permittivity of low-resistivity Si (solid line) sig- 
nificantly departs from the dielectric permittivity of "dielectric" Si in the region 
around the important dimensional parameter of the problem, the characteristic fre- 
quency c/(2a) ~ 10 14 — 10 15 rad/s. That is why the low-resistivity sample cannot 
be described at low frequencies by the static dielectric permittivity of Si equal to 
e Si (Q) = 11.67. To describe it, the term 0(T)/l, like in Eq. ([53]), should be added 
to e Sl (0). Note that in this case (3(T)/l > 1 at the first Matsubara frequencies with 
I = 1, 2, 3, . . . and, thus, this quantity cannot be considered as a small parameter. 
At the same time, for a high-resistivity sample the inclusion of the dc conduc- 
tivity would lead to deviations from the dashed line in Fig. 3 only at frequencies 
£ < 10 s rad/s, which are much less than the characteristic frequency. This compari- 
son permits to make a conclusion in what experimental situations the conductivity 
properties of semiconductors at low frequencies should be taken into account and 
when they should be omitted as being not related to dispersion forces. 

The future expe rime nts on the modification of semiconductor charge carrier 
density by laser lightP^ will bring a more clear understanding of this problem on 
the connection between the low- frequency material properties and dispersion forces. 
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7. DOES SPATIAL DISPERSION LEAD TO AN IMPORTANT 
IMPACT ON THERMAL CASIMIR FORCE? 

The presented above new analytic results on the low-temperature behavior of the 
Casimir free energy, pressure and entropy between dielectrics or between dielectric 
and metal are based on the conventional Lifshitz theory which describes dielectric 
materials by means of the frequency dependent dielectric permittivity. In fact, the 
assumption that the material of the plates possesses only the frequency dispersion 
means that the components of electric displacement are connected with the compo- 
nents of electric field by the relation 

D k {r,uj) = E kl {r,u)Ei{r,uj). (101) 

This equation is central in all different derivations of the Lifshitz formula (see, e.g., 
Refs. H [3 H5H313 [Ml HOI [HI)- The effects of spatial nonlocality (spatial dispersion) 
are in fact essential only at shortest separations between the plates comparable with 
atomic dimensions and also for metals at sufficiently large separations (typically of 
about 2-3 fim) in the frequency region of the anomalous skin effect. The Casimir 
force in the latter region was described by the Lifshitz theory reformulated in terms 
of the Leontovich impedancePS 

Recently the spat ial dispersion c ame to the attention in co nnecti on with the 
thermal Casimir force.^ESESEZISSlin particular, it was claimecP ^^ that for real 
metals at any separation the account of spatial dispersion leads to practically the 
same result ([6]) for the reflection coefficients at zero frequency as was obtained earlier 
using the Drude model dielectric function (|5|). This conclusio n, if it is correct, not 
only returns us to the contradiction with experiment! 1 1 1 but also casts doubts 
on all results obtained by means of the conventional Lifshitz theory accounting for 
only the frequency dispersion. It is natural when the spatial dispersion contributes 
a small fraction of a percent as it is generally believed in numerous applications of 
the Lifshitz theory. It is, however, quite another mat ter when the account of the 
spatial dispersion results in some "dramatic effects" p^i.e., in several hundred times 
larger thermal correction than is obtained in the local case. Below we demonstrate 
that the conclusions of Refs. [64] -[68] are in fact not reliable because they use the 
Lifshitz theory of dispersion forces outside of its application range!^ 

To find the electromagnetic modes associated with an empty gap between the 
plates, Refs. I64l-l68luse the standard continuity boundary conditions, 

E\t — E 2 t, Bi n = B 2n , Di n = D2n, B u — B 2 t, (102) 

which are commonly applied in the derivation of the Lifshitz formula for spatially 
local nonmagnetic materials. Here B is the magnetic induction, n is the normal to 
the boundary directed inside the medium, the subscripts n, t refer to the normal 
and tangential components, respectively, the subscript 1 refers to the vacuum and 
subscript 2 to the plate material. In Refs. [M]- [55] the spatial dispersion is described 
by the longitudinal and transverse dielectric permittivities depending on the wave 
vector and frequency: Ski = w). However, as is shown below, in the theory 
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of the Casimir effect both the boundary conditions (|102[) and permittivities £ki = 
Ski(q,uj) are inapplicable. 

We start from the boundary conditions and recall the set of Maxwell equations 
in a metal describing the Casimir effect, 

1 dB 1 3D 

rotE + — — = 0, div£> = 0, rotB — = divB = 0. (103) 

c ot c at 

Equations (|103p do not contain any external, i.e., independent on E, D and _B, 
current or charge densities. The definition of the electric displacement is 

3D dE 

-W^-at* 4 **' (104) 

where the volume current i is induced by E and B and takes into account the 
conduction electrons. 

In electrodynamics with spatial dispersion the electric field and magnetic induc- 
tion are fi nite at the boundary surfaces, whereas the electric displacement can tend 
to infinity!^ Then, integrating Eqs. (|103[) over the thickness of the boundary layer 
as is done in Ref. [5H1 we reproduce the first two c onditions in Eq. (|102[) and arrive 
at the modified third and fourth conditions jSZESl 

An 

E u = E 2t , B ln = B 2n , D 2n -D ln = 4ttct, [n x (B 2 — B±)] = — j, (105) 

c 

where the induced surface charge and current densities are given by 



±J div[n x [D x n]}dl, j = i- J 



2 dD 



, ........ . ........ , . n -dl. (106) 

Note that the boundary conditions (| 1 05|) , (I106|) are obtained from the macroscopic 
Maxwell equations for physical fields. They should not be mixed with the boundary 
conditions arising in perturbative theories and for the fictitious fields (see below). 

In linear electrodynamics for a medium with time-independent properties with- 
out spatial dispersion the material equation connecting the electric displacement 
and electric field takes the form 

D fc (r,t)= f iu(r,t-lf)Ei(r,lf)dtf. (107) 

According to this equation, the electric displacement at a point r and moment t is 
determined by the electric field at the same point r at different moments t' < t (the 
spatial dispersion is absent but the temporal may be present). It is easily seen that 
the substitution of Eq. (|107p in Eq. (|106p leads to a = 0, j= and, as a result, the 
boundary conditions (|105[) coincide with the standard continuity conditions p02[) . It 
is unjustified, however, to use conditions (|102[) in the presence of spatial dispersion. 
In Refs. [571 [5nia few examples are presented illustrating that in this case neither a 
nor j is equal to zero. 

We now turn to a discussion of the use of dielectric permittivity £ki(q, w) in the 
theory of the Casimir effect with account of spatial dispersion. In the presence of 



5, 2008 2:30 WSPC/INSTRUCTION FILE textDr 



Thermal Casimir force between dielectrics 29 

only frequency dispersion, it is possible to perform the Fourier transformation of 
the fields 

/oo />oo 
E(r,uj)e- lult duj, D(r,t) = D(r,w)e- iut dw, (108) 
- oo J — oo 

in Eq. (]107[) and arrive at Eq. (|101p where 

pOO 

e H (r,w) = / e kl (r,T)e^ T dT (109) 



is the frequency-dependent dielectric permittivity and r = t — t'. In fact Eqs. (|101[) . 
(|108p and (|109p are used in par allel with the bound ary conditions (|102p in all 
derivations of the Lifshitz formulaPE B 6 ! 37 ! 38 ! 69 ! 70 ^ 

If the material of the plates is characterized not only by temporal but also spatial 
dispersion, Eq. (|107p is generalized to 



D k (r,t) = f dt' f dr'e k i(r,r',t-t')Ei{r' 



t'). (110) 



If the material medium were uniform in space (i.e., all points were equivalent), the 
kernel £ would not depend on r and r' separately, as in Eq. (fTT0|) . but on the 
difference R = r — r'. In this case, by performing the Fourier transformation, 

E(r,t) = f duj I dqE(q,Lj)e l(qr - UJt \ 

J —OO J 

D(r,t)= f duj f dqD(q,u J )e^ qr ''^ t \ (111) 

J —OO J 

and substituting it in Eq. (|110[) , one could introduce the dielectric permittivities 

e M (q,u) = £ 'dr J dR£ kl (R,r) e -^ R ^\ (112) 

as Refs. 64-68 do, and rearrange Eq. (IllOp to the form 

D k (q,u) = e k i{q t u)Ei(q,u). (113) 

In the Casimir effect, however, the material medium is not uniform due to the 
presence of a macroscopic gap between the two plates (half spaces). Because of this, 
the assumption that the kernel £ depends only on R and r is wrong. As a result, 
it is not possible to introduce the dielectric permittivity e k i{q,uj) depending on 
the wave vector and frequency. In fact, for systems with spatial dispersion in the 
presence of boundaries the kernel £ depends not only on R and r but also on the 
distance from the boundary!^ In this complicated situati on t he following approx- 
imate phenomenological approach is sometimes applicable. 12 -'- 1 For electromagnetic 
waves with a wavelength A the kernel e(r, r', r) in Eq. (|110[) differs essentially from 
zero only in a certain vicinity of the point r with characteristic dimensions I <C A 
(for nonmetallic condensed media I is of the order of the lattice constant). Then it is 
reasonable to assume that £ is a function of R—r-r', except for a layer of thickness 
I adjacent to the boundary surface. If one is mostly interested in bulk phenomena 
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and neglects the role and influence of a subsurface layer, the quantity Ski{q, ui) may 
be employed as a reasonable approximation. 

This approximate phenomenological approach is widely applied in the theoretical 
investigation of the anomalous skin effect (see, e.g., Ref. l9~Tj) . Note that in Ref. [91] 
some kind of fictitious infinite system was introduced and electromagnetic fields in 
this system are discontinuous on the boundary surface. This discontinuity should 
not be confused with the discontinuity of physical fields of a real system in the 
presence of spatial dispersion described by Eqs. (|105p and (|106|) . (There is also 
another approach to the description of the anomalous skin effect in polycrystals 
using the generalizations of the local Leontovich impedance which takes into account 
the shape of Fermi surface^!) . The frequency and wave vector dependent dielectric 
permittivity in the presence of boundaries is also approximately applied in the 
theory of radiative heat transfei^or in the study of electromagnetic interaction of 
molecules with metal surfaces!^ In all these applications the boundary effects are 
usually taken into account by the boundary conditions (|105j) supplemented by so 
called "additional boundary conditions" . 

It is unlikely, however, that the approximate phenomenological approach using 
such quantity as Eki(q,uj) in the presence of boundaries would be applicable in the 
theory of the Casimir force where the boundary effects on the zero-point electromag- 
netic oscillations are of prime importance. It is notable also that this approach faces 
serious theoretical difficulties including the violation of the law of conservation of 
energyP^Il It is then not surprising that the application of this approach in Refs. [Ml 
1551 results in the contribution to the Casimir free energy fr om the transverse electric 
mode which is in contradiction with experiment EQI32I33] 

One more shortcoming of Refs. 64 -68; is that they substitute the dielectric per- 
mittivity £fc;(<7,cj), depending on both wave vector and frequency, into the conven- 
tional Lifshitz formula derived in the presence of only temporal dispersion. In the 
famous review papeJ^it has been noticed, however, that with the inclusion of spa- 
tial dispersion the free energy of a fluctuating field takes the form T = Tl + AT, 
where Tl is given by the conventional Lifshitz expression derived in a spatially 
local case and written in terms of the Fresnel reflection coefficients, and AT is an 
additional term which can be expressed in terms of the thermal Green's function of 
the electromagnetic field and polarization operator. Review^Sl calls as not reliable 
the results of, e.g., Ref. 97 obtained by the substitution of dielectric permittivity 
£fcj(<7,cij), taking account of spatial dispersion, into the conventional Lifshitz for- 
mula. It can be true that the Lifshitz formula written in terms of general reflection 
coefficients is applicable in both spatially local and nonlocal situations. However, 
as far as the exact reflection coefficients in a spatially nonlocal case are unknown, 
the use of some approximate phenomenological models, elaborated in literature for 
applications different than the Casimir effect, may lead to incorrect results for AT 
and create inconsistencies with experiment. 
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To conclude this section, the results of Refs. [64}|68| on the influence of spatial 
nonlocality on the Casimir interaction are shown to be not reliable. Although at 
present there is no fundamental theory of the thermal Casimir force incorporating 
spatial dispersion, there is no reason to expect that it can play any significant role 
in the frequency region of infrared optics (experimental separations) or normal skin 
effect (i.e., at separations between plates greater than 4-5 /mi). 

8. CONCLUSIONS AND DISCUSSION 

In the foregoing we have presented the derivation of analytic asymptotic expres- 
sions for the free energy, pressure and entropy of the Casimir interaction between 
two dielectric plates and between metal and dielectric plates at low and high temper- 
atures. It was shown that the low-temperature behavior of the Casimir interaction 
between dielectrics and between dielectric and metal is determined by the static 
dielectric permittivities of nonpolar dielectrics. The obtained results were shown to 
be in agreement with thermodynamics when the static dielectric permittivities of 
dielectrics are finite. In particular, the entropy of the Casimir interaction goes to 
zero when the temperature vanishes, i.e., the Nernst heat theorem is satisfied. This 
demonstrates the consistency of the original Lifshitz's approach to the van der Waals 
forces between dielectrics which disregards the small conductivity of dielectrics at 
constant current. 

The second important result shown above is that the inclusion of the dc con- 
ductivity of dielectrics into the model of dielectric response leads not to some small 
corrections to the characteristics of the Casimir interaction, as one could expect, 
but makes the Lifshitz theory inconsistent with thermodynamics leading to the vi- 
olation of the Nernst heat theorem. This reveals that real material properties at 
very low, quasistatic frequencies are in fact not related to the phenomenon of the 
van der Waals and Casimir forces which is actually determined by sufficiently high 
frequencies. In this case the zero-frequency contribution to the Casimir force should 
be understood not literally but as analytic continuation to zero of the material 
physical behavior in the region around the characteristic frequency. 

The presented results provide a basis for the calculation of the van der Waals 
and Casimir forces between real materials. Such calculations are much needed for 
numerous applications of the Casimir force discussed in the Introduction, in par- 
ticular for the applications in nanotechnology and for constraining predictions of 
fundamental physical theories beyond the Standard Model. Bearing in mind that 
semiconductors are the main constituent materials in nanotechnological devices, it 
is a subject of high priority to understand the Casimir effect with semiconductor 
boundaries. In this connection we have discussed new experimental results and the- 
oretical ideas on the Casimir interaction between a metal sphere and semiconductor 
plate. It was stressed that by changing the charge carrier density in the semicon- 
ductor it is possible to bring it in different intermediate states between metallic 
and dielectric. In this case the problem arises when the conductivity properties of 
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semiconductor are not related to dispersion forces and when they are becoming rel- 
evant. A criterion for the resolution of this problem was formulated based on the 
relation between the typical frequency at which the dc conductivity properties come 
into play and the characteristic frequency of the Casimir effect. In fact the thermal 
Casimir interaction between semiconductors remains an open question and much 
work should be done both in experiment and theory to gain a better insight into 
this subject. 

The last major problem discussed in this review is whether or not the spatial 
dispersion influences essentially the thermal Casimir force between real materials. In 
the present literature there is no agreement on this subject. We adduced arguments 
in favor of the statement that in the region of experimental separations the influence 
of the spatial dispersion on the Casimir force is negligible small. The statements 
on the opposite, contained in the literature, were shown to be not reliable because 
they are obtained by the application of the Lifshitz theory outside of its application 
range. At the same time it was ascertained that at the moment there is no consistent 
fundamental theory of the van der Waals and Casimir forces taking the spatial 
dispersion into account. This is the problem to solve in the foreseeable future. 
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